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Abstract. The majority of VO-compatible spectra handling applications operates only with 
a few spectra entirely downloaded from single or several SSAP servers. We try to identify 
the scientific cases which could immediately benefit from future SSAP applications de- 
signed for GRID deployment. Their key feature is the sophisticated spectra pre-selection 
and preprocessing done on distributed servers using the intelligent agent summarising the 
results and performing final high-level processing or analysis. 
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1. Introduction 

There is a wealth of calibrated astronomi- 
cal spectra accessible in current data archives, 
however most of them are not suitable for di- 
rect physical analysis. Being pre-reduced by 
some automatic pipeline or individually by 
manual reduction, they are mostly stored in 
archives in form of FITS file or ASCII table 
as the relation of intensity in arbitrary numbers 
(instrumental counts, digital numbers etc.) and 
wavelength (or frequency). 

The scientific analysis of such spectra re- 
quires further processing by the variety of dif- 
ferent methods. In certain studies a huge num- 
ber of spectra has to be collected from different 
servers (e.g. in different spectral regions ) and 
transformed into common units. 

The clear definition of metadata descrip- 
tion and easy unit conversion as well as trans- 
parent lookup and download of spectra in 
Virtual Observatory has a great potential to be- 



come a new advanced way of astronomical re- 
search. The connection of powerful distributed 
computing available by GRID, Web Services 
and VO protocols can establish an innovative 
research environment allowing the data mining 
of huge datasets. 



2. Current status of VO tools for 
spectroscopy 

Although there is a number of various tools 
written for the analysis of spectra of astro- 
nomical objects (e.g IRAF packages splot 
and spectool, Starlink DIPSO, MIDAS 
XALICE as well as numerous single-purpose C, 
FORTRAN and IDL routines), the number of 
modern tools understanding the VO protocols 
(especially spectra access protocol - SSAP) is 
still very low. We are giving a brief overview 
of their capabilities and status of development 
(in July 2008). 
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2.1. VOSpec 

http : //esavo . esa . int/vospec/ 

- Developed at ESA 

- Very simple (mainly for building SED) 

- Polynomial and Gaussian fits only 

- Blackbody fit 

- No RV measurement built-in 

- It can be called directly from VizieR 

- Can work with linelists through Simple 
Line Access Protocol (SLAP) 

- Theoretical VO (TVO) supported (syn- 
thetic spectra: Kurucz, disks) 

- Rapid development 

- Integration of some methods in users own 
program possible 

- JAVA applet 

- Contains PLASTIC VO-interoperability 
layer 

- De-reddening of extragalactic objects 
built-in 

- Support of dimensional equations for units 
(DIMEQ, SCALEQ) 

- Rather complicated view of users data - 
needs to create SSA wrapper and prepare 
VOTable even for one simple FITS spec- 
trum before viewing 

2.2. SpecView 

http : //www . stsci . edu/resources/ 
software_hardware/specview 

- It is available as JAVA applet or stand-alone 
application 

- Developed and supported by STScI 

- Understands a number of formats from 
HST instruments and most NASA/ESA 
satellites 

- Easy work with local data (binary or ASCII 
tables and ID FITS) 

- Well tailored for practical spectral analysis 

- Powerful linelists (many included) 

- Fitting of models by x 2 — even user mod- 
els may be constructed 

- Built-in number of standard stars spectra 
and models for various physical conditions 

- Whole Kurucz library of spectra available 
for immediate display. 



- Only simple polynomials for visual fit 
available 

- De-reddening, CLOUDY models built-in 

- Proper error propagation of fitted values 

- Support of dimensional equations 

- Does not support PLASTIC 

2.3. SPLAT-VO 

http : //star-www . dur . ac . uk/"pdraper/ 
splat/splat-vo/ 

- Supported by JCMT after closing Starlink 

- Most advanced for stellar astronomy (fit- 
ting by freehand drawing, INTEP) 

- RV measurement - both Gaussian fit and 
mirroring 

- Custom line list for individual spectra 
available 

- Wavelet analysis 

- Full featured data and editor and spread- 
sheet 

- Publication quality output, powerful plot- 
ting options, annotations 

- Supported PLASTIC 

- Reads ID FITS image files (e.g. rebinned 
IRAF multispec files) 

- Asynchronous SSAP queries 

2.4. Period04 

http : //www . univie . ac . at/tops/PeriodQ4/ 

- Its predecessor was widely used Period98 

- Period04 is rewritten in Java (and partly C) 

- Needs formatted text files 

- supports VOTable files 

- Several period finding methods built-in 

- Can handle even the regular period shifts 

- PLASTIC layer under preparation 

2.5. FROG 

http : //www . j ach . hawai i . edu/sof tware/ 
starlink/ 

- Part of Starlink Java package 

- Included in last JCMT Starlink release 
(Hokulei) 

- VO-protocols built-in 

- Almost same capabilities as Period04 
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- Not developed after closing Starlink 

- Last version since 2004 

- The period searching engine available re- 
motely as a SOAP service 

- Web service for Fourier transform 

- Collaborates with TOPCAT, but not 
PLASTIC 

- Easy but powerful for doing period analy- 
sis in VO environment 



- (De)convolution of instrumental profile 

- Application of physical broadening func- 
tions (rotation, limb darkening) 

- Shift in radial velocity, application of he- 
liocentric correction computed on server 

- Merging of individual echelle orders in a 
single spectrum 



2.6. SDDS Spectrum Services 

http : //voservices . net/spectrum 

- Advanced access and visualization of 
Sloan spectra 

- Search of similar spectra using PCA 
method 

- Several continuum and line fitting proce- 
dures 

- Programmable interface available but in 
.Net or Mono framework 

- Collaboration with other SSA servers not 
supported - Only SDSS data at JHU 

3. Spectra post-processing 

The current VO applications using SSA proto- 
col are very simple providing only the URL of 
the given dataset, so the client gets a whole (of- 
ten quite large) file. The further processing is 
then done fully by the client who has to down- 
load all the relevant datasets and keep them in 
the memory or local storage. The real strength 
of the VO technology lies, however, in trans- 
ferring part of the client's work to the server, 
which usually runs on a powerful machine 
with fast connection to the data archives. The 
typical example of commonly required post- 
processing of fully reduced spectra (at least in 
stellar astronomy) are given below: 

- Cutout services (selection of only certain 
spectral lines or regions giving the wave- 
length range) 

- Projection of multidimensional datasets (in 
3D spectroscopy) 

- Rectification of continuum 

- Rebinning to given, usually equidistant 
grid of wavelengths (constant AA or A In A) 



Although most of these post-processing 
methods can be implemented in a straightfor- 
ward way, the rectification of continua may be 
a difficult art or black magic. Especially two 
problems are most common and most difficult 
to solve easily. 

Normalising echelle spectra 
A really challenging problem (which had not 
yet been solved fully) is the normalisation 
of echelle spectra - especially with wide ab- 
sorption lines or complicated profiles (P Cyg, 
Be stars). Some wide lines may span sev- 
eral echelle orders and so the individual or- 
ders have to be precisely unblazed and then 
merged. This problem is very complex, de- 
pending on construction of individual spec- 
trograph, the observation strategy etc. The 
explanation of pro blem and partly the solu - 
tion is given by ISkoda & Hensbereel d2003[). 
ISkoda & Slechtal d2004. iPrugniel & Soubiranl 



(120011) or lErspamer & North! (2002) 



X 1 fitting procedures 

The more advanced way of normalising is the 
least squares fitting based on some theoreti- 
cal model. It is very popular in radioastron- 
omy and X-ray astronomy (XSPEC) rather for 
estimation of the global SED fit parameters 
(usually only as combination of several sim- 
ple power law functions). Quite useful in op- 
tical spectroscopy is the usage of appropri- 
ate synthetic spectrum to find the real contin- 
uum in data. The implementation of such a 
procedure in VO tools could call the Theory 
VO (TVO) servers to get correct models (e.g. 
Kurucz) and, if necessary, to improve the mod- 
els by recomputing on the GRID in a iterative 
way. Unfortunately, the models of many inter- 
esting objects are very speculative. 
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Fig. 1. Dynamical spectrum of H/3 line profile 
variability of 59 Cyg. Residuals from average 
profile of 38 spectra are binned to 20 phase 
bins corresponding to period 28. 192 days. Two 
expanded cycles shown for clarity. Ind ividual 
profile s are overplotted above. After iMaintj 
(120031) 



4. Advanced spectral analysis 

As the goal of VO is to make easy and comfort- 
able the physical analysis of a huge number of 
fully reduced (and post-processed) spectra in 
the environment of VO client (or web portal), 
the common recipes of spectral analysis have 
to be implemented as VO-compatible. While 
the simplest ones may be done by client, the 
full power of VO may be exploited only by 
their implementation on server side (and best 
deployed on GRID due to their inherent in- 
dependencies). Here are the examples of most 
commonly used techniques of spectra analysis 
in the stellar astronomy: 



4.1. Dynamical spectrum 

It is sometimes called the gray representation 
or trailed spectrum. The basic idea is to find the 
small time-dependent deviations of individual 
line profiles from some average. 



First the average of many high dispersion 
high SNR spectra (with removal of outliers) 
is prepared (called template spectrum). Then 
each individual spectrum in time series is either 
divided by the template (quotient spectrum) or 
the template is subtracted from it (the differ- 
ential spectrum). The group of similar result- 
ing intensities is given the same colour or level 
of gray. See Fi g. [1] More examp l es may 
be fou nd e.g. in | de Jong et all dl99 9), Maintz 
(120031) or lUvtterhoevenl(l200l . 

The result is drawn in 2D image where on 
horizontal axis is a wavelength in the line pro- 
file or corresponding RV relative to laboratory 
wavelength, on the vertical axis the time of 
middle of observation (in HJD) or the circu- 
lar phase when the data are folded with certain 
period. 

4.2. Measurement of radial velocity and 
higher moments of line profile 

The one of the important information received 
from spectrum is the radial velocity (RV) of 
the object. From its changes the binarity can 
be revealed, or the possession of extrasolar 
planet. RV is sometimes presented as zero-th 
moment of the spectral line profile. The higher 
moments of the line profile are important as 
well. The first moment is equivalent width. The 
combination of higher moments of line pro- 
file is a one of the possible ways of determi- 
nation of non radial pulsa tion modes - num- 
bers I, m dAertset al.ll 19921) . The RV is usually 
measured from single lines (by finding their 
centre using fits or mirroring) or from part of 
spectrum by cross-correlation techniques using 
the template spectra (observed or synthetic) of 
similar object. 



4.3. Measurement of equivalent width 

The Equivalent width (EW) of the spectral line 
gives the information about the number of ab- 
sorbing or emitting atoms of given element. 
The measurement of equivalent width (EW) of 
a given line needs the determination of area 
bordered by line profile and the continuum. 
This area is transformed to the rectangular line 
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with the depth of 1 .0. The width of this rectan- 
gle is called equivalent width. Emission lines 
by definition have negative EW. The changes in 
EW during time may bring about the informa- 
tion about the dynamic evolution of the target 
and may be subjected to period analysis. 

Due to its dependency on continuum place- 
ment, a problem of automatic rectification 
arises again. It is extremely difficult on echelle 
spectra and spectra of late type stars. If not 
rectified properly, the extreme ly shallow lines 
may give large error in EW (lVrancken et al.1 
119971) . 

EW are often used in abundance analy- 
sis. A good check of correct normalisation is 
the comparison of EW of spectral lines of 
the same element in different spectral regions 
dHensberge et alJ2000l) . 

Some attempts of auto matic rectification 
have been proposed eg. by IZhao et al.l (2006) 
but in general another kind of information 
(synthetic spectrum or template from other in- 
strument) is required. 

The measurement of EW of a number of 
lines in thousands of spectra may be accom- 
plished as a VO service and run in parallel on 
a GRID. 



4.4. Bisector analysis 

It is a method describing quantitatively the tiny 
asymmetry or subtle changes in line profiles. It 
is easily done by marking the middle of hor- 
izontal cuts of profile in different line depths. 
The line connecting such points (called bisec- 
tor) is then zoomed in horizontal (wavelength 
or RV) direction. See Fig. [2] 

The characteristic shape of bisector gives 
the information about turbulence fields (e.g. 
convection) in stellar photosphere, charac- 
teri sed by the value of micro-turbulent veloc- 
ity dGravfll982ll2005l) or about other processes 
causing the tiny profile asymmetry. It has been 
used su ccessfully for searc hing of extrasolar 
planets dPovich et al.l 1200 ll) or in asteroseis- 
mology. Requires high resolution (echelle) and 
high SNR normalised spectra. 
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Fig. 2. The construction of bisector of line pro- 
file (left pan el) and the zo omed bisector (right 
panel). From lGravl d!982h . 



4.5. Period Analysis 

One of the most popular methods used in as- 
tronomy is the period analysis. Its aim is to 
find the hidden periods of variability of given 
object. Sometimes this period can be identi- 
fied with some physical mechanism (e.g. or- 
bital period of binaries, rotational modulation 
or pulsations). Wide range of objects show the 
multi-periodicity on various time scales (e.g. 
binary with pulsating components). Ones the 
suspected period is found, the data may be 
folded accordingly, plotted in circular phase 
corresponding this period. Very helpful is a 
interactive capability of showing data folded 
while selecting different peaks at the peri- 
odogram. 

The current availability of VO tools is lim- 
ited. The nice FROG (part of Starlink new 
Java-based development) is abandoned and so 
the only hope is the Period04 by R Lenz, who 
is working on the PLASTIC compatible data 
exchange of VO tables (Lenz, 2008, private 
communication). 

Periodogmm of line profiles. It is a kind of 
gray representation with period on vertical axis 
and the zoomed line profile on horizontal. 
See Fig. [3] Works well for anal ysis of multi- 
perio dical non-radial pulsations dde Jong et alj 
1999). 
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Fig. 3. Example of 2D periodogram of Silll 
4552 A line profile variability of A Sco. The 
darker gray level marks the higher value of 
power spectr um in given wavelen gth for given 
period. After U vtterhoevenl (12004 . 



4.6. Doppler imaging 



It was introduced by IVogt & Penrodl {l983a) 
as a method allowing the surface mapping of 
stellar spots. First test were d one on stars of 
RS CV n type and on f Oph dVogt & Penrodl 
1983b). Works well on rapid rotators and needs 
a high resolution spectra with very high SNR 
(300-500). Tho whole rotational period should 
be covered well, better several times. When all 
the requirements are met, the map of surface 
features (spots, nodes of non radial pulsations) 
is obtained with very high accuracy. See the 
left panel of Fig. [4] 



4. 7. Zeeman Doppler Imaging 

Quite complicated processing of spectra is re- 
quired for study of stellar magnetic fields. The 
estimation of magnetic field from polarimetry 
using the Zeeman phenomena involves the pro- 
cessing of long series of homogeneous spec- 
tra to be accomplished in parallel with ex- 
treme precision and requires again the in- 
formation from synthetic models (simulation 
of Stokes parameters on simulated magnetic 
sta rs) The nice exam ple is the model of II Peg 
by IStrassmeierl d2008l) . See the right panel of 
Fig. a 




Fig. 4. Zeeman Dopple r Imaging of II Peg 
After lStrassmeierl |2P08) 



4.8. Spectra disentangling 

This method allows to separate the spectra 
of individual stars in binary or multiple sys- 
tems and simultaneously to find orbital pa- 
rameters of the system, even in case of heavy 
blending of lines. It supposes the changes in 
line profile are caused only by combination of 
Doppler shifted components (no intrinsic vari- 
ability of star). The best solution of orbital pa- 
rameters and disentangled line profiles of in- 
dividual stellar components are found by least 
square global minimisation. The method also 
enables to remove the telluric lines with great 
precision. The good orbital coverage and the 
estimate of orbital parameters is required. Two 
approaches exist: 

Wavelength space disentangling 
developed bvlBagnuolo & Giesl(ll99ll) and im- 
proved by ISimon & SturrrJ (119941) . It needs a 
large memory to store sparse matrices, requires 
large computing power. It is more straightfor- 
ward to understand the results and sources of 
errors. 

Fourier space disentang l ing 
introduced by iHadraval d 1 995L 1 19971) in pro- 
gram KOREL. Another program available 
today (still based on KOR EL ideas) is 
FDBINARY dlliiic et alJbOOH) . They work in 
Fourier space, and transform the wavelengths 
into In A. They solve a small amount of linear 
equations, so they are memory savvy and can 
be run on even small computer. The method, 
however, requires perfect continuum fit (dif- 
ficult to achieve at merged echelle spectra). 
On unreliabl y rectified spectr a strange artifacts 
may appear (Iliiic et al. 2004). 
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4.9. Deconvolution of Iodine lines 

One of the new approaches to high precision 
radial velocity measurement is the usage of 
Iodine vapour cells in front of the spectro- 
graph. The idea is to use the huge amount of 
sharp Iodine lines as a ruler for marking exact 
position of investigated stellar line thus elimi- 
nating the instrumental shifts. 

However, the Iodine lines have to be 
first removed from spectra to see the stel- 
lar ones. The process is not so straight- 
forward and requires the models and tem- 
plates from high resolution Fourier Transform 
Spectrographs as well as exposures without 
the cell. Finally the global optimization (e.g. 
by maximum entropy method) i s required. See 
e.g. Martinez Fiorenzand (120061) . 

4.10. Classification of stellar spectra 

It is a one of the methods where the full power 
of the Virtual observatory infrastructure could 
be exploited to work on massive data vol- 
umes. The goal is to assign to every spectrum 
a short string of letters and numbers repre- 
senting the physical characteristics of the star 
(surface temperature, size, mass, peculiar be- 
haviour etc.). 

It used to be a time consuming method 
done by visual comparison of many spectra 
from photographic plates (e.g. MKK classifica- 
tion). Today it has been done finding the min- 
imal differences between a grid of template 
(or even synthetic) spectra and the examined 
one. The techniques of artificial intelligence 
(e.g. neural networks, genetic algorithms, self- 
organizing maps, cluster analysis techniques) 
may be very efficient, but the global optimiza- 
tion using the x 1 minimization seems to be 
more reliable with possibility of iterative con- 
trol. An example of automatic classi fication 
engine for white dwarfs is described bv l Winter! 
d2006l) . 

As it can be seen, the algorithm can run 
in parallel and uses workflows, so it may be 
nicely implemented as Grid application of VO 
services. 



5. Spectral analysis in Grid 
supported VO services 

5.1. Workflows and pipelines 

Many astronomical legacy applications are de- 
signed to run as a single task on local PC, the 
exchange of results with other tools is done 
over the local files and a processing is driven 
by a set of parameters. Some packages have to 
be run as separate tasks in given order as some 
task are dependent on a information written to 
the processed file by the preceding tasks. 

So the work has to be done in a strictly 
serial manner — step by step. Thus the idea 
of workflows is quite natural for the most of 
scientist's work but there is a strong reluc- 
tance of astronomers (driven by a strong con- 
servatism and habits) to use modern work- 
flow builders/executors in the data reduction or 
spectra analysis. The certain level of workflow 
management is accepted for distributed com- 
puting of models on clusters — usually using 
batch queue scheduler scripts either on GRID 
or on institutional network of PC (e.g. Condor, 
Beowulf). 



5.2. Massive parallel spectra reduction 

The processing of spectra coming from current 
sophisticated spectrographs is a quite computer 
demanding task and thus may be the good 
candidate for processing on GRID. Especially 
the rigorously correct reduction of echelle 
spectra requirin g the non-linear Le a st S quares 
optimization (Piskunov & Valentil l2002l) or 
even the estimate of bac kground light model 
dChurchil& Allenlll995l) for the hundreds of 
echelle exposures would benefit from GRID 
as it can be run in parallel for all frames 
using the same calibrations (flat, bias) and 
templates (order position, aperture estimates). 
Some pipeplines are already written using the 
concept of workflows (e.g. Gasgano, Reflex for 
ESO UVES and FORS spectrograph, respec- 
tively). 
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5.3. Advantages of VO services 

The key VO feature, that allows the novel ap- 
proach to the everyday scientific work and thus 
promises new discoveries, is the interoperabil- 
ity of all VO-compatible tools and common 
format for describing metadata (VOTable). 

It allows (in principle) the most of current 
work of an astronomer to be done in one GUI 
which can provide an transparent access to VO 
tools accomplishing all of the most tedious and 
routine work: 

The VO infrastructure enables transparent 
search, download and on-the-fly conversion of 
data, unified presentation of different relations 
from multi-dimensional datasets, the compli- 
cated reduction and analysis tasks to be done 
in background or deployed on GRID, the re- 
sults may be stored directly in databases or 
staged on virtual network home space (like 
MySpace). The most sophisticated data pro- 
cessing may be implemented in the form of 
web services (WS). 

Although most of the methods described 
may be (and have been) run in local envi- 
ronment using the legacy astronomical appli- 
cation, their implementation as VO services 
would benefit from the VO infrastructure and 
might allow the analysis of a huge amount of 
data easily. The key issues in VO approach are 
the 

- Unified data format (VO-Table, semantics 
of variables) 

- Quick search for available data in a global 
manner using SSAP and registries 

- Transparent data conversion, homogeniza- 
tion, rescaling 

- Powerful presentation with remote data 
and theoretical results (TVO) 

- Staging disks with virtual home (VOSpace, 
MySpace) 

Combining the VO infrastructure power and 
the easy and transparent high performance 
computing on GRID will allow the advanced 
analysis of large spectral surveys feasible in a 
reasonable time. It may force even the kind of 
serendipitous research — e.g. you may click 
on star in the image of globular cluster to see 
its dynamic spectra folded with estimated pe- 
riod. 



6. Killer spectral applications for VO 

Here we give some scientific cases that may 
be the killer applications forcing the stellar as- 
tronomers to use the VO tools suggested above 
despite their current reluctance and hesitation: 

- Use VO to find all stars with emission in 
given line, find the time when it was in 
emission and plot the time evolution of its 
EW. 

- Use VO to get 1000 or more spectra of the 
given object, cut out regions around given 
lines, plot the lines, make a gray dynamic 
spectrum folded in time 

- The same, but fold by period clicked on in- 
teractive periodogram 

- Get the unknown lines identification of 
piece of spectra from theoretical observa- 
tory (Simple line access protocol) having 
the line selection limited by pre-estimated 
temperature (using Saha equation) 

- Create light and radial velocity curve of a 
binary star for given period (estimated by 
other VO tool running in parallel and ex- 
changing data over PLASTIC) 

- Fit the grid of models (r eff , log g) to the ob- 
served spectrum for many stars (e.g. from 
given cluster) 

7. Conclusions 

Astronomical spectroscopy uses a wide range 
of techniques with different level of complexity 
to achieve its final goal — to estimate the most 
precise and reliable information about celestial 
objects. The large part of spectroscopic anal- 
ysis today has been accomplished by several 
independent non VO-compatible legacy pack- 
ages, where each works with different local 
files in its own data format. Analysis of large 
number of spectra is thus very tedious work re- 
quiring good data bookkeeping. 

Accomplishing the analysis in VO infras- 
tructure may benefit from automatic aggre- 
gation of distributed archive resources (e.g. 
the multispectral research), seamless on-the-fly 
data conversion, common interoperability of 
all tools (using PLASTIC protocol) and pow- 
erful graphical visualisation of measured and 
derived quantities. 
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The currently available VO clients support- 
ing the SSA protocol can provide with only ba- 
sic functions in a interactive environment. The 
more advanced work with spectra (e.g. auto- 
matic preprocessing or advanced analysis) is 
not supported in VO at all. The most of the 
legacy applications should be turned into the 
VO server side services and the conversion of 
legacy scripts and recipes into Workflows will 
allow their easy deployment on GRIDs. 

Modern techniques of astronomical analy- 
sis require the considerable amount of comput- 
ing power and very offten the iterative com- 
parison with theoretical models is indispens- 
able. The properly designed VO service get- 
ting both observed and synthetic spectra trans- 
parently from VO resources and processing the 
number of objects in parallel on GRID might 
become a killer application turning the atten- 
tion of wide astronomical community to the 
Virtual observatory as a viable and innovative 
way of modern astronomical research. 

By introduction of modern VO-aware tools 
into the astronomical spectral analysis a re- 
markable increase of effectiveness of astro- 
nomical research can be achieved. 
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